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摘要 :气候 和 土地 利用 变化 是 影响 水 资源 变化 最 直接 的 因素 。 应 用 SWAT 模型 对 干旱 半 干 旱 区 小 流域 宁夏 清水 河流 域 径流 进 
行 多 情景 模拟 预测 ,以 历史 气候 要 素 变化 趋势 和 CA-Markov 模型 分 别 设置 未 来 气候 和 土地 利用 变化 情景 ,以 决定 系数 尺 和 
Nash-Sutcliffe 模型 效率 系数 Ens ( Nash-Sutcliffe efficiency coefficient) 来 衡量 模拟 值 与 实测 值 之 间 的 拟 合 度 ,并 评价 模型 在 清水 
河流 域 的 适用 性 。 结 果 表 明 , 昔 府 湾 站 在 校准 期 和 验证 期 的 RR 分 别 为 0.80 和 0.71,Ens 分 别 为 0.77 和 0.69 ,泉眼 山 站 在 校准 期 
和 验证 期 的 尼 分 别 为 0.66 和 0.63,Ens 分 别 为 0.62 和 0.56, 表 明 构 建 的 SWAT 模型 可 以 用 于 清水 河流 域 的 径流 模拟 。 对 未 来 
气候 和 土地 利用 变化 情景 下 径流 的 模拟 结果 显示 ,径流 变化 主要 由 降水 变化 主导 ,降水 减少 和 气温 升 高 的 综合 作用 对 流域 径流 
影响 最 为 显著 ; 由 于 耕地 和 建设 用 地 的 增加 ,未 来 3 种 土地 利用 情景 下 流域 径流 量 将 均 会 呈现 明显 增加 变化 。 与 2010 年 
相 比 ,到 2020 年 ,自然 增长 情景 流域 径流 将 增加 17.04% ,林地 保护 情景 径流 将 增加 14.44% ,规划 情景 径流 将 增加 13.98% ;综合 
降水 .气温 和 土地 利用 的 结合 变化 情景 显示 ,未 来 流域 径流 将 会 有 不 同 程度 的 下 降 ,规划 情景 和 气候 变化 的 结合 4 0 
降 最 为 明显 ,而 有 意 增 大 林地 和 加 强生 态 保护 的 林地 保护 情景 对 减缓 流域 径流 下 降 具 有 一 定 作 用 。 在 气候 变化 的 大 背景 下 

据 水 资源 利用 管理 目标 ,可 通过 调整 流域 管理 措施 ,特别 是 土地 利用 变化 和 改善 区 域 小 气候 来 减缓 气候 变化 对 流域 水 资源 
面 效果 ,以 此 来 改善 流域 径流 和 生态 环境 状况 。 
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Abstract: Climate and land use changes are very important factors that directly influence water resources. To predict runoff 
under climate and land use changes in Ningxia Qingshui River Basin, which is a typical basin located in the arid and semi- 
arid area of Yellow River Valley, the distributed hydrologic SWAT model was calibrated and validated using the coefficient 
of determination ( R’*) and Nash-Sutcliffe coefficient ( Ens) at the Hanfuwan and Quanyanshan hydrological stations. Several 
future climate and land use change scenarios for the year 2020 were set using the climate historical trend and CA-Markov 


model. The calibration and validation results for Hanfuwan station showed that R’ was 0.80 and 0.71, and ENS was 0.77 and 
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0.69, respectively, while the results for Quanyanshan station showed that R* was 0.66 and 0.63, and ENS was 0.62 and 
0.56, respectively. The results suggested that the SWAT model could be used to predict the runoff under climate and land 
use changes in Qingshui River basin. Further, our results showed that runoff was dominated by precipitation, and the result 
under the scenario of reduced precipitation combined with increased temperature had the largest impacts. Runoff can be 
expected increase in the future land use scenarios because of the increasing farmland and built-up area. Under natural 
growth, planned development, and woodland/ ecological protection scenario, the runoff in 2020 would increase by 17.04%, 
14.44% , and 13.98% respectively compared with 2010. Lastly, the runoff change in 2020 was predicted by SWAT model 
under the combine of land use scenarios and climate change scenarios. The results showed that the runoff would decrease in 
the future under all combinations of scenarios. Among them, the runoff decreased perceptibly under the combination of 
planned development scenario and climate changes, and the woodland/ ecological protection measure is helpful to slow down 
the rate of runoff decline under all climate change scenarios in the future. Therefore, in the context of climate change, 
ecologically prudent constructions and land-use managements are necessary to cope with the severe situation of runoff change 


in Qingshui River basin. 
Key Words: climate change; land use change; Ningxia Qingshui River Basin; SWAT model; runoff prediction 


环境 变化 下 的 水 文 研究 是 当今 全 球 变 化 研究 的 热点 之 一 1。 在 环境 变化 中 ,气候 变化 和 土地 利用 变化 
是 影响 水 资源 变化 最 直接 的 两 个 因素 。 气 候 变 化 直接 改变 大 气 降 水 的 空间 分 布 和 时 间 变 异 特性 中, 土地 利 
用 变化 则 能 直接 引起 地 表 蒜 散 ( 发 ) 变 化 从 而 使 水 资源 发 生 改 变 。 由 于 分 布 式 水 文 模型 既 考 虑 了 流域 的 空间 
异 质 性 ,也 能 较 好 地 对 流域 水 文物 理 过 程 进行 刻画 ,因而 被 广泛 运用 于 两 者 变化 对 水 资源 影响 的 研究 当 
中 1* 引 其 中 ,SWAT( 水 土 评 价 工具 ,Soil and water assessment tool) 模 型 已 被 大 量 研究 证 明 其 在 各 种 地 理 和 气 
候 区 域 均 具 有 较 好 的 适用 性 !“7” 。 尤 其 在 水 资源 贫乏 的 干旱 半 干 旱 区 ,SWAT 模型 已 经 成 为 水 文 研究 最 重要 
的 工具 之 一 |。 

我 国 干旱 半 干 旱 区 面积 广阔 ,但 自然 条 件 恶 劣 ,水 资源 量 少 质 差 , 且 空 间 分 布 严重 不 均 。 人 类 的 社会 经 济 
活动 基本 集中 在 黄河 流域 和 内 陆 河流 域 片区 ,使 得 流域 周边 人 口 密度 大 ,人 均 占 有 量 严重 不 足 ,水 资源 已 经 成 
为 干旱 半 干 旱 区 社会 经 济 发 展 的 主要 制约 因素 。 对 该 区 域内 流域 的 水 文 变 化 特征 ,尤其 是 在 气候 和 土地 利用 
变化 背景 下 的 未 来 径流 变化 特征 进行 深入 研究 ,并 将 研究 结果 应 用 于 流域 水 资源 规划 与 管理 中 ,这 对 构建 健 
康 稳 定 的 干旱 半 干 旱 区 流域 生态 系统 极为 重要 。 在 目前 已 有 的 基于 SWAT 模型 的 气候 和 土地 利用 变化 下 水 
文 变化 的 研究 中 ,未 来 气候 变化 情景 的 常见 方法 有 IPCC 法 "和 基于 历史 变化 趋势 变化 法 “中 两 种 。 由 
于 IPCC 预测 结果 的 不 确定 性 及 较 低 的 分 辨 紊 等 局 限 性 ,因此 基于 气候 变化 趋势 设置 变化 情景 的 方法 在 相关 
研究 中 得 到 更 为 广泛 的 应 用 。 但 在 大 多 气候 变化 情景 研究 中 ,人 研究 者 通常 将 流域 作为 一 个 整体 来 对 待 ,并 未 
考虑 流域 的 空间 异 质 性 ,使 得 研究 结果 的 实用 价值 受到 影响 。 而 在 未 来 土地 利用 变化 情景 的 设置 中 ,也 多 以 
土地 利用 现状 为 基础 ,简单 设置 不 同 土地 利用 类 型 面积 的 增 减 变 化 , 较 少 结合 研究 区 域 的 实际 和 发 展 规划 。 

宁夏 清水 河 是 宁夏 境内 最 大 的 一 条 黄河 支流 ,其 中 多 年 平均 地 表 水 可 利用 量 仅 为 0.7 亿 mx ,但 流域 内 
近年 取水 总 量 超 过 4 亿 mi ,是 黄河 上 游 典 型 的 缺 水 型 流域 。 近 年 来 由 于 地 下 水 位 严重 下 降 , 矿 化 度 上 升 ， 
流域 生态 安全 面临 巨大 挑战 。 本 文 尝试 基于 SWAT 模型 构建 流域 的 分 布 式 水 文 模 型 ,通过 气候 要 素 变化 趋势 
分 区 设置 气候 变化 情景 ,并 利用 CA-Markov 模型 预测 未 来 土地 利用 变化 情景 ,以 此 分 析 两 者 对 流域 产 流 的 影 
响 ,并 对 流域 未 来 河流 径流 进行 模拟 预测 。 人 研究 结果 可 为 宁夏 清水 河流 域 水 资源 保护 利用 提供 参考 ,并 为 黄 
河流 域 乃 至 干旱 半 干 旱 区 域 的 水 资源 管理 .生态 环境 保护 .政策 制定 以 及 合理 规划 提供 决策 依据 与 科学 指导 。 


1 研究 区 概况 


宁夏 清水 河 发 源 于 六 盘山 东 入 固原 市 原 州 区 开 城 镇 黑 刺 沟 脑 ,向 北 流 经 固原 市 的 原 州 区 \ 西 吉 县 ,中 卫 市 


http ://www.ecologica.cn 


201703.00059v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1254 生 态 学 报 37 卷 


的 海原 县 . 沙 坡 头 区 .中 宁县 和 吴忠 市 的 同心 县 ` 红 寺 堡 
区 等 区 县 ,于 泉眼 山西 侧 注 入 黄河 ,全 长 303 km。 在 宁 
夏 境 内 流域 地 理 位 置 介 于 35$"53' 一 37"28'N ,105"12' 一 
106*38 中 之 间 ,面积 13511 km? ,是 宁夏 境内 最 大 的 黄 
河 一 级 支流 (图 1)。 流 域 地 势 南 高 北 低 , 地 形 地 貌 复 
森 , 以 黄土 丘陵 为 主 , 中 上 游 为 洪 积 一 冲积 平原 区 ,中 游 
西 侧 为 黄土 丘陵 、 盆 . 霄 区 ,中 游 东 侧 为 黄土 丘陵 区 ,中 
下 游 为 河谷 平原 。 属 温带 半 干 旱 气 候 区 ,多 年 平均 降水 
量 为 335 mm, 时 空 分 布 不 均匀 ,主要 集中 在 7 一 9 月 ( 占 
全 年 降水 量 的 70% 以 上 ), 自 下 游 至 上 游 降 水 在 200 一 
600 mm 之 间 ; 年 均 气 温 6.2%C ; 多 年 平均 蒸发 量 1000 一 
1400 mm。 


37°00N 


36°00’ 


2 研究 区 SWAT 模型 构建 


2.1 模型 输入 数据 106°00'E 
SWAT 模型 的 输入 数据 包括 DEM 数据 (90 mx90 图 1 宁夏 清水 河流 域 位 置 示意 图 
mm, 来源 于 国际 科学 数据 服务 平台 ) i 土壤 数据 土地 利 Fig.l Location of Ningxia Qingshui River basin 


用 数据 、 气 象 数 据 和 水 文 数据 。 土 壤 数 据 包括 土壤 空间 
分 布 数据 和 土壤 物理 属性 数据 ,物理 属性 数据 中 土壤 机 械 组 成 数据 需 使 用 MATLAB 中 采用 3 次 样 条 插值 法 
将 国际 制 转换 为 美国 制 ' 引 。 土 地 利用 数据 以 流域 2000 年 和 2010 年 两 期 TM 影像 (来 源 于 国际 科学 数据 服务 
平台 ) 解 译 获得 ,结果 经 检验 Kappa 系数 分 别 达 到 0.81.0.78。 气 象 数据 为 流域 内 部 及 周边 6 个 气象 站 点 
1980 一 2010 年 的 日 观测 数据 (来 源 于 中 国 气象 数据 网 ) 。 水 文 数据 为 清水 河流 域 中 游 韩 府 湾 水 文 站 和 下 游 录 
眼 山 水 文 站 2000 一 2011 年 的 月 径流 实测 数据 。 
2.2 ”模型 校准 及 验证 

由 于 ArcSWAT 模型 自 带 的 模型 参数 敏感 性 分 析 、 模 型 校准 和 验证 工具 的 功能 相对 较 弱 "9 ,因此 选取 
SWAT-CUP 对 清水 河流 域 SWAT 模型 进行 模型 校准 和 验证 。 将 2000 年 作为 预 热 期 ,2001 一 2005 年 为 校准 
期 ,2006 一 2011 年 为 验证 期 。 遵 循 流 域 自 上 而 下 的 校准 原则 , 先 对 韩 府 湾 水 文 站 以 上 区 段 ,后 对 韩 府 湾 至 泉 
眼 山 间 区 上段 进行 参数 敏感 性 分 析 和 校准 。 
2.3 ”模型 适用 性 评价 

选取 决定 系数 RR 和 Nash-Sutcliffe 模型 效率 系数 Ens 来 衡量 模拟 值 与 实测 值 之 间 的 拟 合 度 中 ,以 此 评 
价 SWAT 模型 在 清水 河流 域 的 适用 性 。 

决定 系数 尼 表 示 模 拟 值 与 实测 值 的 吻合 程度 ,表达 式 ! 9 为; 


2 


之 (Q,; - Q%) (0,; - Qs) 


尺 = 
[> (0 = 0 ”> CO = i) | 全 
式 中 ,0, ;为 实测 流量 值 ,0, ,为 模拟 流量 值 , 0, 为 实测 流量 平均 值 , 0, 为 模拟 流量 平均 值 。 有 R 取 值 为 
0 一 1, 数 值 越 小 表示 吻合 程度 越 低 , 一 般 尽 达 到 0.50 以 上 模拟 结果 可 以 接受 ,0.70 以 上 模拟 比较 准确 ,等 于 1 
表明 非常 吻合 [2 。 
Ens 能 判断 水 文 模型 的 模拟 效果 ,表达 式 0 为 : 
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> (Q,, 
Ens=1-— 
> (0,; - 0%.)? 


式 中 ,0, ;为 实测 流量 值 ,0, ;为 模拟 流量 值 , 0, 为 实测 流量 平均 值 。Ens 在 0 一 1 之 间 变 动 , 越 接近 于 1 模拟 
精度 越 高 ,一 般 在 0.50 一 0.65 之 间 为 满意 ,0.65 一 0.75 之 间 为 好 ,0.75 一 1.00 之 间 为 很 好 。 


3 ”情景 设置 


3.1 气候 变化 情景 

根据 宁夏 气候 区 划 '”3 ,结合 流域 地 理 特征 ,将 清水 河流 域 划分 为 北部 .中 部 和 南部 3 个 分 区 ,北部 包括 中 
卫 .中 宁 两 个 气象 站 ,中 部 为 同心 站 ,南部 包括 海原 . 西 吉 和 固原 3 个 气象 站 。 基 于 1980 一 2010 年 3 个 分 区 气 
象 站 的 降水 和 气温 变化 趋势 ,以 2010 年 气象 数据 ( C0 ) 为 基础 ,以 10a 为 研究 尺度 ,设置 2020 年 气候 变化 情景 
如 下 : 

情景 1(C1) 只 降水 发 后 改变 。 北 部 汛期 降水 增加 0.87 mm , 非 汛 期 减少 2.57 mm; 中 部 汛期 下 降 11.84 
mm , 非 汛期 降水 下 降 4.68 mm; 南 部 汛期 降水 下 降 10.05 mm , 非 汛 期 减少 2.34 mm。 

情景 2(C2) 只 气温 发 生 改 变 。 北 部 最 高 温 上 升 0.49% ,最 低温 上 升 0.59% ;中 部 最 高 温 上 升 0.38%C ,最 
低温 上 升 0.39% ;南部 最 高 温 上 升 0.59C ,最 低温 上 升 0.71%C 。 

情景 3(C3) ”降水 和 气温 同时 发 生 改变 的 组 合 情 景 。 
3.2 土地 利用 变化 情景 

根据 流域 历史 土地 利用 变化 情况 ,结合 未 来 发 展 规划 ,构建 流域 2010 一 2020 年 的 3 种 土地 利用 发 展 情 
景 ,利用 CA-Markov 模型 对 流域 2020 年 的 土地 利用 进行 模拟 ,并 将 2010 年 实际 土地 利用 设置 为 基础 情景 
(S0) ,以 此 来 分 析 未 来 土地 利用 变化 对 径流 的 影响 (图 2) 。 


语言 计 


0 20 40km 
国王 林地 加 水 域 国 国 未 利用 地 流域 边界 
加 草地 大 建设 用 地 [耕地 


2 清水 河流 域 2020 年 土地 利用 类 型 空间 分 布 图 
Fig.2 Distribution of land use types in Qingshui River basin of 2020 
S0:; 基础 情景 Base scenario; S1 :自然 增长 情景 Natural growth scenario; S2: 规划 情景 Planning scenario; S3 :林地 保护 情景 Woodland/ecological 


protection scenario 


情景 1(S1) 为 自然 增长 情景 ,土地 利用 变化 按 原 有 趋势 持续 演变 。 
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情景 2(S2 ) 为 规划 情景 ,在 2010 年 土地 利用 的 基础 上 ,参考 宁夏 2006 一 2020 年 土地 利用 总 体 规划 ,调整 
流域 2020 年 各 土地 利用 类 型 的 面积 。 

情景 3(S3) 为 林地 保护 情景 ,在 规划 情景 的 基础 上 ,严格 限制 其 他 土地 利用 类 型 侵占 林地 和 水 域 ,对 坡度 
大 于 25" 耕 地 全 部 进行 退耕 还 林 的 举措 。 
3.3 气候 和 土地 利用 结合 变化 情景 

气候 变化 和 土地 利用 变化 协同 作用 于 径流 变化 ,因此 设置 两 者 结合 变化 情景 ,如 表 1 所 示 。 


表 1 流域 气候 和 土地 利用 结合 变化 情景 


Table 1 Scenarios of climate and land use combination 


土地 利用 变化 情景 气候 变化 情景 Climate change scenario 
Land use change scenario Cl C2 C3 
Sl SI1+C1 S1+C2 S1+C3 
S2 S2+C1 S2+C2 S2+C3 
S3 S3+C1 S3+C2 S3+C3 


ma S1 :自然 增长 情景 Natural growth scenario; S2: 规划 情景 Planning scenario; S3 :林地 保护 情景 Woodland/ecological protection scenario; C1 :降水 
变化 情景 Precipitation change scenario; C2; 气 温 变 化 情景 Temperature change scenario; C3 : 降水 与 气温 结合 变化 情景 Precipitation and temperature 


combination scenario 


4 结果 与 分 析 


4.1 SWAT 模型 校准 及 验证 

选择 和 径流 有 关 的 14 个 参数 进行 敏感 性 分 析 ,综合 清水 河上 段 和 下 段 的 参数 分 析 结 果 可 知 对 清水 河流 
域 径流 产生 最 敏感 的 参数 有 CN2( SCS 模型 径流 曲线 参数 ,SCS runoff curve number for moisture condition II)、 
CH_N2( 主 河道 曼 宁 系数 ,Manning's "n" value for the main channel) .CANMX (最 大 覆盖 度 参 数 , Maximum 
canopy storage) 和 ALPHA_BNK( 河 岸 蕃 水 基 流 a 系数 ,Base flow alpha factor for bank storage ) 等 ,其 余 参 数 影 
响 相 对 较 小 ( 表 2)。CN2 的 变化 直接 影响 径流 产生 ,CH_N2 的 大 小 对 流域 汇流 时 间 影 响 显 著 , 从 而 进一步 影 
响 着 水 文 过 程 线形 状 的 变化 。 径 流产 生 对 CANMX 的 变化 敏感 说 明 清 水 河流 域 径流 产生 受 植被 覆盖 变化 影 
天 响 较 大 ,径流 变化 对 ALPHA_BNK 的 变化 敏感 说 明 清 水 河流 域 河岸 的 蓄 水 能 力 在 对 径流 产生 影响 也 较 大 。 
(= 在 SWAT-CUP 工具 中 ,模型 参数 校准 与 参数 敏感 性 分 析 同 时 进行 , 韩 府 湾 及 泉眼 山水 文 站 的 校准 结果 如 
也 。 表 3 所 示 。 韩 府 湾 水 文 站 月 径流 量 的 RR 为 0.80，Ens 为 0.77, 泉 眼 山水 文 站 的 民 和 Ens 分 别 为 0.66 和 0.62， 
表明 模型 在 清水 河流 域 月 径流 的 校准 期 结果 较为 满意 。 

为 检验 模型 在 清水 河流 域 的 适用 性 ,对 校准 后 的 模型 进行 验证 。 由 表 3 可 知 韩 府 湾 水 文 站 在 验证 期 的 
RR 和 Ens 分 别 达 到 0.71 和 0.69 ,模拟 效果 较 好 。 泉 眼 山 水 文 站 月 径流 在 验证 期 的 尼 为 0.63，Pms 为 0.56, 模 
拟 结果 仍 在 可 接受 范围 。 

结果 表明 ,构建 的 SWAT 模型 可 以 用 于 清水 河流 域 的 河流 径流 模拟 , 且 在 韩 府 湾 以 上 区 段 的 适用 性 要 好 
于 其 以 下 区 段 。 流 域 上 段 和 下 段 产生 适用 性 差异 是 由 模型 校准 过 程 所 决定 ,SWAT 模型 校准 属于 分 段 校准 ， 
在 上 段 进行 校准 后 再 对 下 段 进 行 校 准 , 上 段 校准 模拟 的 径流 结果 直接 用 于 下 段 河 流 径 流 的 模拟 ,从 而 导致 模 
型 在 下 段 的 模拟 精度 下 降 。 

4.2 情景 模拟 
4.2.1 气候 变化 下 流域 径流 模拟 

利用 SWAT 模型 对 清水 河流 域 出 口 录 眼 山 站 2020 年 河流 径流 量 对 气候 变化 的 响应 进行 模拟 ,结果 见 表 
4。 仅 降水 发 生变 化 时 ,径流 变化 响应 较为 敏感 ,2020 年 年 均 径流 量 下 降 到 1.197 m:/s, 相 比 2010 年 产 流下 降 
3.82% ;而 气温 升 高 变化 对 径流 变化 影响 相对 较 弱 ,径流 量 下 降 了 0.022 ms, 变 化 率 仅 为 -1.82%; 当 降 水 和 
气温 同时 发 生变 化 时 , 即 流域 降水 减少 ,温度 也 升 高 的 情况 下 ,径流 变化 响应 非常 明显 ,径流 变化 率 达 到 
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-5.13% ,年 均 径 流量 仅 为 1.180 m/s。 


表 2 模型 参数 敏感 性 分 析 结 果 


Table 2 Sensitivity analysis of the parameters 


韩 府 湾 Hanfuwan 泉眼 山 Quanyanshan 
排序 ”参数 定义 排序 ”参数 定义 
Order Parameter Parameter definition Order Parameter Parameter definition 
1 CN2 SCS 径流 曲线 参数 1 CN2 SCS 径流 曲线 参数 
2 CH_K2 沟 道 有 效 水 导电 率 2 CH_N2 主 河道 曼 宁 系数 
3 CH_N2 主 河道 曼 宁 系 数 3 CANMX 最 大 和 覆盖 度 参数 
4 CANMX 最 大 覆盖 度 参数 4 ALPHA_BNK 河岸 蓄 水 基 流 a 系数 
5 ALPHA_BNK 河岸 蓄 水 基 流 a 系数 5 SOL_AWC 土壤 可 利用 水 量 参 数 
6 SOL_K 饱和 水 力 传导 参数 6 SLSUBBSN 平均 坡 长 参数 
7 SOL_AWC 土壤 可 利用 水 量 参数 7 SOL_Z 土壤 深度 参数 
8 ESCO 土壤 蒸发 补偿 参数 8 GWQMN 浅 层 地 下 水 径流 参数 
9 REVAPMN 浅 层 地 下 水 再 蒸发 参数 9 REVAPMN 浅 层 地 下 水 再 蒸发 参数 
10 ALPHA_BF 基 流 消退 参数 10 ESCO 土壤 蒸发 补偿 参数 
11 GWQMN 浅 层 地 下 水 径流 参数 11 CH_K2 沟 道 有 效 水 导电 率 
12 SOL_Z 土壤 深度 参数 12 ALPHA_BF 基 流 消退 参数 
13 SLSUBBSN 平均 坡 长 参数 13 SOL_K 饱和 水 力 传 导 参 数 
14 BLAI 最 大 潜在 叶 面 积 指 数 参数 14 BLAI 最 大 潜在 叶 面 积 指数 参数 


于 


CN2: SCS runoff curve number for moisture condition II; CH_K2. Effective hydraulic conductivity of main channel; CH_N2: Manning's "n" value for 
the main channel; CANMX : Maximum canopy storage; ALPHA_BNK: Base flow alpha factor for bank storage; SOL_K: Saturated hydraulic conductivity ; 
SOL_AWC: Available water capacity of the soil layer; ESCO: Soil evaporation compensation coefficient; REVAPMN : Threshold water level in shallow 
aduifer for revap; ALPHA_BF: Baseflow recession constant; CWQMN : Threshold water level in shallow aquifer for base flow; SOL_Z: Depth from soil 
surface to bottom of layer; SLSUBBSN : Average slope length; BLAI: Potential maximum leaf area index for the plant 


表 3 清水 河流 域 SWAT 模型 月 径流 校准 和 验证 结果 
Table 3 Calibration and validation of SWAT model for monthly runoff in Qingshui River basin 
月 均值 Monthly mean value/ (m3/s) 


id cg validation 、 实测 值 模拟 值 R ns 
Observation Prediction 

校准 期 Calibration period(2001 一 2005) 韩 府 湾 0.81 0.68 0.80 0.77 

泉眼 山 2.43 2.57 0.66 0.62 

验证 期 Validation period(2006 一 2011) 韩 府 湾 0.41 0.37 0.71 0.69 

泉眼 山 1.33 1.60 0.63 0.56 


Ens: Nash-Sutcliffe 模型 效率 系数 Nash-Sutcliffe efficiency coefficient 


表 4 不 同 气 候 变 化 情境 下 流域 径流 变化 


Table 4 Changes of runoff under climate scenarios in Qingshui River basin 


气候 变化 情景 年 均 径 流量 /(m’/s) 径流 变化 率 /% 气候 变化 情景 年 均 径 流量 /(m3/s) 径流 变化 率 /% 
Climate change scenario ”Annual runoff value Change rate of runoff ||Climate change scenario ”Annual runoff value Change rate of runoff 
C0 1.244 0.00 Cl 1.197 —3.82 
C2 1;222 —1.82 C3 1.180 —5.13 


C0: 气候 基础 情景 Climate base scenario; C1: 降水 变化 情景 Precipitation change scenario; C02. 气温 变化 情景 Temperature change scenario; C3 : 


降水 与 气温 结合 变化 情景 Precipitation and temperature combination scenario 


总 体 而 言 ,流域 未 来 径流 变化 的 过 程 中 ,降水 与 气温 结合 变化 对 径流 的 影响 最 大 ,其 中 降水 下 降 引 起 的 径 
流 减少 占 主导 作用 ,这 一 结果 与 西北 地 区 的 其 他 相关 研究 结果 相同 ”1 。 降 水 变化 直接 影响 流域 的 水 资源 总 
量 , 其 变化 对 径流 变化 具有 确定 作用 ,气温 变化 影响 流域 的 燕 散 ( 发 ) 变 化 , 随 着 温度 升 高 ,流域 内 的 车 发 增 
加 ,在 降雨 不 变 情况 下 ,径流 减少 ” ,而 在 当 蒸发 水 量变 化 不 大 时 ,降水 对 径流 影响 更 显著 '*。 
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4.2.2 土地 利用 情景 下 流域 径流 模拟 

未 来 3 种 土地 利用 情景 下 径流 变化 模拟 结果 如 表 $ 所 示 , 到 2020 年 3 种 土地 利用 情景 下 的 径流 量 较 
2010 年 均 有 明显 上 升 。 变 化 最 大 的 是 自然 增长 情景 ( 31) ,径流 增加 了 17.04% , 主要 原因 是 流域 耕地 和 建设 
用 地 有 较 大 程度 的 增加 ( 表 6) ,清水 河流 域 河谷 平原 区 耕地 增加 受 限 ,耕地 的 增加 转变 为 山地 .丘陵 区 的 坡 耕 
地 的 增加 ,造成 径流 的 显著 增 大 六 ,此 外 ,建设 用 地 的 增长 增 大 了 流域 不 透水 面 , 有 助 于 径流 增加 。 规 划 情 景 
(S2) 的 径流 量 相 对 S1 情景 有 所 下 降 , 但 对 于 2010 年 仍 有 13.98% 的 变化 率 ,可 能 是 S2 情景 下 耕地 增加 虽然 
明显 ,但 建设 用 地 .未 利用 地 和 草地 相对 S1 情景 有 明显 减少 所 致 ; 林地 保护 情景 (S3 ) 相 对 于 S2 情景 , 林 


=. 


地 水域 等 有 所 增加 ,耕地 和 建设 用 地 有 所 减少 ,导致 流域 径流 稍 有 增加 "3 。 


表 5 不 同 土地 利用 情景 下 流域 径流 变化 


Table S Changes of runoff under land use scenarios in Qingshui River basin 


土地 利用 情景 


土地 利用 情景 


年 均 径 流量 / (m/s) 径流 变化 率 /% 年 均 径 流量 / (m/s) 径流 变化 率 /% 
Land use Land use 
Annual runoff value Change rate of runoff Annual runoff value Change rate of runoff 
change scenario change scenario 
S0 1.244 0.00 S1 1.456 17.04 
S2 1.418 13.98 S3 1.424 14.44 


S0: 基 础 情景 Base scenario; S1 :自然 增长 情景 Natural growth scenario; S2 : 规划 情景 Planning scenario; S3 : 林地 保护 情景 Woodland/ecological 


protection scenario 


表 6 清水 河流 域 不 同 土地 利用 情景 下 土地 利用 类 型 面积 
Table6 Structure of land use types of Qingshui River basin in 2020 


S0 Sl S2 S3 
土地 利用 类 型 
Land use type 面积 /km? 比例 /% 面积 /km ”比例 /% ”面积 /km ”比例 /% 面积 /km? 比例 /% 
Area Percentage Area Percentage Area Percentage Area Percentage 

林地 Woodland 409.10 3.11 409.11 3.11 451.45 3.43 467.80 3.55 
草地 Grassland 6641.77 50.44 5479.01 41.61 5412.08 41.10 5412.08 41.10 
水 域 Water 165.18 1.23 164.49 1.25 164.35 1.25 166.51 1.26 
建设 用 地 Building land 402.62 3.06 586.87 4.46 449.16 3.41 447.90 3.40 
未 利用 地 Unused land 119.19 0.91 52.19 0.40 6.42 0.05 6.42 0.05 
耕地 Farmland 5429.34 41.23 6475.54 49.18 6683.74 50.76 6666.49 50.63 


S0: 基 础 情景 Base scenario; S1 :自然 增长 情景 Natural growth scenario; S2 : 规划 情景 Planning scenario; S3 : 林地 保护 情景 Woodland/ecological 


protection scenario 


结果 显示 坡 耕 地 的 增加 和 生态 用 地 的 增加 均 能 对 流域 径流 产生 积极 的 影响 ,但 是 在 坡 耕 地 增加 的 同时 流 
域 水 土 流失 和 土壤 侵蚀 问题 也 可 能 进一步 加 剧 ,因此 进行 植树 造林 退耕 还 林 还 草 等 增加 生态 用 地 的 方法 才 
是 未 来 流域 土地 利用 的 转变 方向 。 

4.2.3 气候 与 土地 利用 结合 变化 下 流域 径流 模拟 

对 气候 和 土地 利用 变化 协同 作用 下 流域 2020 年 的 径流 进行 模拟 ,结果 显示 ( 表 7) 在 降水 减少 (Cl ) 情况 
下 ,未 来 3 种 土地 利用 情景 的 流域 径流 均 呈 下 降 变化 ,变化 率 介 于 -0.67% 一 -3.41%; 只 考虑 气温 变化 时 
(C2) ,3 种 土地 利用 情景 的 流域 径流 均 有 较 大 程度 的 增加 (12.31% 一 15.38% ) ;降水 与 气温 结合 作用 情景 
(C3) ,流域 径流 则 呈现 不 同 程度 的 下 降 , 其 下 降 变化 率 介 于 -1.62% 一 -4.24% 之 间 , 均 略 大 于 C1 情景 。 在 不 
同 气候 变化 情景 条 件 下 ,3 种 土地 利用 情景 的 径流 变化 也 均 表 现 出 S1> S3> S2 ,结果 与 径流 量 对 土地 利用 变 
化 响应 的 结果 一 致 。 

综合 而 言 ,在 模拟 的 9 种 组 合 类 型 中 ,最 有 可 能 发 生 的 变化 情景 是 土地 利用 、 降 水 和 气温 的 综合 变化 情 
景 , 即 SI1+C3、S2+C3 和 S3+C3 情景 ,其 中 S2+C3 情景 径流 下 降 变化 最 大 ,变化 率 为 -4.24% ，S3+C3 次 之 
(-3.87% ) ，S1+C3 情景 径流 下 降 变化 最 小 (-1.62% ) 。 说 明 相 对 于 S2+C3 情景 ,有 意 增 大 林地 和 加 强生 态 保 
护 的 S3+C2 情景 在 未 来 对 减缓 流域 径流 下 降 具 有 一 定 作用 。 
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表 7 气候 与 土地 利用 结合 变化 下 流域 径流 变化 


Table7 Changes of runoff under combination scenarios in Qingshui River basin 


土地 利用 变化 气候 变化 Climate change scenario 
Land use change scenario C0 C1 C2 C3 
S0 1.244 1.197 1.222 1.180 
年 均 径 流量 / (m/s) Sl 1.456 1.236 1.435 1.224 
Annual runoff value S2 1.418 1.202 1.397 1.191 
S3 1.424 1.207 1.403 1.196 
土地 利用 变化 气候 变化 Climate change scenario 
Land use change scenario C0 C1 C2 C3 
S0 0.00 -3.82 -1.82 -5.13 
径流 变化 率 /% Sl 17.04 -0.67 15.38 -1.62 
Change rate of runoff S2 13.98 -3.41 12.31 一 4.24 
S3 14.44 -3.01 12.75 -3.87 


S0; 基 础 情景 Base scenario;S1: 自 然 增长 情景 Natural growth scenario; S2: 规划 情景 Planning scenario; S3 : 林地 保护 情景 Woodland/ecological 


protection scenario; C0: 气候 基础 情景 Climate base scenario; Cl: 降水 变化 情景 Precipitation change scenario; 02; 气温 变化 情景 Temperature change 


scenario; C3 : 降水 与 气温 结合 变化 情景 Precipitation and temperature combination scenario 


在 气候 变化 的 背景 下 ,人 为 规划 干扰 下 的 土地 利用 变化 情景 的 径流 下 降 最 为 明显 ,而 有 意 增 大 林地 和 加 
强生 态 保护 的 规划 对 减缓 流域 径流 下 降 具 有 一 定 作 用 。 土 地 利用 方面 要 从 生态 的 角度 增加 流域 径流 ,因此 要 
在 减少 和 改善 丘陵 山区 坡 耕 地 的 同时 ,加快 林地 ,草地 等 生态 用 地 的 建设 。 而 在 改善 小 气候 上 仍 依托 于 土地 
利用 变化 ,地 表 林 地 、 草 地、 水 域 等 均 在 一 定 程度 上 能 够 改善 区 域 小 气候 ,改善 土壤 微 环境 ,降低 地 表 温 度 , 减 
少 水 分 的 直接 燕 发 ,对 地 表 径 流产 生 具 有 积极 作用 。 


S 结论 与 讨论 


构建 宁夏 清水 河流 域 SWAT 模型 ,对 模型 的 适用 性 进行 评价 ;建立 几 种 未 来 气候 变化 和 土地 利用 变化 情 
景 ,利用 SWAT 模型 对 流域 2020 年 径流 变化 进行 模拟 ,得 到 以 下 几 点 结论 ; 

构建 的 SWAT 模型 可 以 适用 于 清水 河流 域 的 河流 径流 模拟 , 且 在 韩 府 湾 以 上 的 上 游 区 段 适用 性 较 好 ,在 
韩 府 湾 以 下 的 中 下 游 区 段 适用 性 有 所 下 降 。 以 后 的 研究 可 通过 继续 细 化 和 深化 收集 数据 来 增加 SWAT 模型 
的 模拟 精度 ,提高 模型 的 适用 性 。 

流域 径流 量 在 构建 的 几 种 气候 变化 情景 中 均 呈 减 小 变化 ,径流 量变 化 由 降雨 变化 主导 ,降水 减少 和 气温 
升 高 的 综合 作用 对 流域 径流 影响 最 为 显著 。 流 域 径流 量 在 未 来 几 种 土地 利用 情景 下 均 将 会 较 大 程度 地 增 大 
变化 ,其 中 自然 增长 情景 下 流域 径流 有 最 大 程度 的 增加 ,其 次 为 林地 保护 情景 和 规划 情景 。 

流域 2020 年 最 有 可 能 发 生 的 3 种 土地 利用 降水 和 和 气温 的 结合 情景 下 的 模拟 结果 显示 流域 径流 均 会 有 
不 同 程度 的 下 降 , 即 未 来 清水 河流 域 河流 径流 极 可 能 发 生 下 降 变化 。 由 于 清水 河流 域 及 宁夏 黄河 流域 的 流域 
生态 系统 处 于 干旱 胁迫 之 中 ,河流 径流 的 减少 变化 使 得 农业 、 生 态 需 水 等 用 水 形势 更 加 严峻 ,也 将 对 生态 环 
境 造成 较 大 影响 。 

在 气候 变化 的 大 背景 下 ,根据 水 资源 利用 管理 目标 ,可 通过 调整 流域 管理 措施 ,特别 是 土地 利用 变化 和 改 
善 区 域 小 气候 等 减缓 气候 变化 对 流域 水 资源 的 负面 效果 。 因 此 ,清水 河流 域 在 未 来 流域 管理 和 规划 过 程 中 ， 
要 从 生态 保护 的 角度 出 发 ,通过 调整 土地 利用 类 型 ,合理 布局 ,从 而 减缓 径流 的 下 降 , 甚 至 增 大 流域 出 水 量 。 
本 研究 也 可 为 黄河 流域 及 干旱 半 干 旱 区 其 他 流域 的 水 资源 管理 .生态 环境 保护 政策 制定 以 及 合理 规划 提供 
理论 参考 。 
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